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We discuss an earlier application of the Ritz variational method for strongly nonlinear 
problems. We clearly prove that the results derived for several extremely simple problems 
of supposedly physical and mathematical interest do not provide any clue on the utility 
of the approach. 

1. Introduction 

In a recent review in this journal, He [1] analyzed several asymptotic methods for 
strongly nonlinear equations. One such approach, the so called Ritz method, consists 
mainly in converting the nonlinear differential equation into a Newton equation of 
motion. Thus, by minimization of the action integral for the Lagrangian function 
one obtains an approximate solution to the nonlinear equation that is expected to be 
optimal from a variational point of view. Obviously, the accuracy of the approximate 
solution depends heavily on the chosen variational ansatz or trial "trajectory" . 

The purpose of this paper is to analyze the Ritz method proposed by He [1] and 
determine if it is already useful for solving actual nonlinear problems. For the sake 
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of clarity in what follows we devote a section to each of the problems discussed. 

2. Anharmonic oscillator 

First, He [1] transforms the equation of motion for the Dufhng oscillator 

u{t) - u{t) + eu{tf = (1) 

into the variational integral 

2 2 4 / ^ ' 

and concludes that "it requires that the potential V{u) = + must be 

positive for alH > 0, so an oscillation about the origin will occur only if eA^ > 2, 
where A is the amplitude of the oscillation". Unfortunately, He [1] does not show 
the trial function from which he draws that conclusion. 

The Duffing oscillator has been widely studied and, consequently, its properties 
are well known [2] . For example, from straightforward inspection of the potential 
V{u) we already know that there is an unstable equilibrium point at u = when 
e < 0. On the other hand, when e > the potential V{u) exhibits a local maximum 
y = at M = and two minima of depth — l/(4e) symmetrically located at il/y^. 
If the initial conditions are such that V{A) > then the oscillation will certainly be 
about the origin; otherwise there will be oscillations about one of the two minima. It 
is clear that by simple inspection of the potential we obtain much more information 
that the one derived by He [1] from the action integral. Therefore, He's application 
of the Ritz method to this model is of no relevance whatsoever. 

3. A chemical reaction 

He [1] proposed the application of the Ritz method to the chemical reaction 

nA^C + D (3) 

If Nj^{t), NB{t), and Nc{t) arc the number of molecules of the species A, B, and C, 
respectively, at time t then He [1] assumed that Na{0) = a, and Nb{0) = Nc{0) = 0. 
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If we call X = NB{t) = Nc{t), then we conclude that NA{t) = a ~ nx, where x is 
known as the extent of reaction [3,4]. The unique rate of reaction can be defined in 
terms of the extent of reaction as v = dx/dt. He [1] further assumed that the rate 
law is given by 

^^Ha^xr (4) 

At this point we stress the fact that this expression is correct only if the chemical 
reaction ([3]) is elementary, otherwise the rate law may be more complicated. Most 
chemical reactions are not elementary and therefore the reaction order and molecu- 
larity do not necessarily agree, as discussed in any book on physical chemistry [3] or 
chemical kinetics [4] . What is more, the order of reaction may not even be a positive 
integer [3,4]. For concreteness here we assume that the rate law ^ is correct. 

He [1] obtained an approximate solution to the differential equation ^ by means 
of the action integral 

1 



' 2 



dx^ 
Hi 



+ k\a - x) 



2n 



dt (5) 



and the variational ansatz 

where is a variational parameter. Notice that x'"°''^{t) satisfies the boundary con- 
ditions aX t — and t oo. He [1] found that the optimal value of the effective 
first-order rate constant r/ was given by 

V = —r- (7) 

Furthermore, He [1] argued that chemists and technologists always want to know the 
half-time — t{x — a/2) (which he called halfway time). According to equations 
([6]) and ([7]) the half-time is given approximately by 

ivar 

According to He [1] the exact reaction extent for n = 2 is 

.f-(n = 2).a(l-^) (9) 
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This result is obviously wrong because it exhibits an unphysical pole dX t — 1/ {ka). 
From this incorrect expression He [1] derived a meaningless negative half-time 

tf/2(n = 2) = -^ (10) 

In order to obtain a reasonable agreement with the variational result ^ He [1] then 
carried out the following wrong calculation 

(11) 

^'^^ ' ka"-^ ka ^ ' 

In this way He [1] managed to obtain two unphysical negative half-times that agreed 
98%. 

Disregarding the mistakes outlined above we may ask ourselves whether the 
approximate variational result may be of any utility to a chemist. Any textbook 
on physical chemistry [3] or chemical kinetics [4] shows that the exact solution to 
equation (jH) is 

^e^act ^ ^ K \ \ ,n^l (12) 

\ [l + fc(n-l)a"-i<]'/("-'^ J 
and that the exact half-time is given by 

"fc(n-l)a"-i ^ 
It is common practice in chemistry to estimate the half-time from experimental data 
in order to determine the order of the reaction. Obviously, an inaccurate expression 
would lead to an inexact order of reaction. 

The variational half-time ([5]) is reasonably accurate for n = 2 because it is exact 
for n = 1. The reason is that the variational ansatz ^ is the exact solution for 
a first-order reaction when rj ~ k. Notice that equation ([7]) leads to such a result 
when n = 1. We can easily verify that the ratio i™2/*i/2'^* increasingly deviates 
from unity as n increases. Therefore, n = 2 (the only case selected by He [4]) is the 
most favorable case if n is restricted to positive integers greater than unity. 

The half-time (or half-life) is a particular case of partial reaction times. We 
may, for example, calculate the time t = ti/4 that has to elapse for the number of 
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A molecules to reduce to a/4 {x = 3a/4). It is not difficult to verify that 

^exact 

ll±- = 2"-i + 1 (14) 

4.exact ^ ' 

^1/2 

From the experimental measure of ti/2 and ti/^ chemists are able to obtain the 
reaction order n. However, if they used He's variational expression ^ they would 
obtain 

j-var 

^ = 2 (15) 

4-var ^ ^ 

4/2 

that is useless for n ^ 1. According to what we have said above it is not surprising 
that this ratio is exact for n — 1. We clearly appreciate that the variational result 
does not provide the kind of information that chemists would like to have because 
it only predicts first-order reactions. 

From the discussion above we conclude that no chemist will resort to the varia- 
tional expressions in the study of chemical reactions. There is no reason whatsoever 
for the use of an unreliable approximate expression when one has a simple exact an- 
alytical one at hand. Besides, we have clearly proved that the variational expressions 
are utterly misleading. 

4. Lambert equation 

He [1] also applied the Ritz method to the Lambert equation 

y"{x) + —y{x) = (1 - n)^-^ (16) 
and arrived at the variational formulation 

J{y) = \ I {-nV'"'y" + kV") dt (17) 

By means of the transformation z = He obtained 

J{z) = ]- [ + k^z^) dt (18) 



2 _ 

that leads to the Euler-Lagrange equation 



z" + k^z = (19) 
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Obviously, the solution to this linear equation is straightforward. 

If we substitute the transformation z — y"^ into equation (jl6p we obtain equation 
(|19p in a more direct way. Therefore, there is no necessity for the variational Ritz 
method. 



5. Soliton solution 

He [1] also studied the KdV equation 

du{x,t) du{x,t) d^u{x,t) 

__ _ 6^.(., t)^— + = 

and looked for its travelling-wave solutions in the frame 



(20) 



u{x,t) ^u{£,), i^x- ct (21) 
The function u(^) satisfies the nonlinear ordinary differential equation 

u"'{0 - cu'iO - QuiOu'iO - (22) 

where the prime indicates differentiation with respect to ^. Then He [1] integrated 
this equation (taking the integration constant arbitrarily equal to zero) and obtained 



u"{0 - cuiO - MO' = 



(23) 



By means of the so called semi-inverse method He [1] obtained the variational 
integral 



J = 



Choosing the trial function 



1 / du 



2 



-u +u 



dt 



u — pcosh ^(g^) 



(24) 



(25) 



where p and q are variational parameters. He [1] obtained p — c/2 and q — \/c/2. 

By substitution of equation (|25l) into equation ([22ll we obtain the same values 
of p and (7 in a more direct way and with less effort. Therefore, there is no need for 
the variational method for the successful treatment of this problem. 
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6. Bifurcation 

He [1] also applied the Ritz method to the most popular Bratu equation 

u"(a;) + Ae"(^)=0, u(0) = u(l) = (26) 

that has been studied by several authors [5] (and references therein). Here we only 
cite those papers that are relevant to present discussion. He [1] derived the action 
integral 

J = Qw'^ - Ae"^ dx (27) 

and proposed the simplest trial function that satisfies the boundary conditions: 

u{x) = Ax{l - x) (28) 

Curiously, He [1] appeared to be unable to obtain an exact analytical solution 
for the integral; however, it is not difficult to show that 

J (A) = ^ - erf ( V^/2) (29) 

We cannot exactly solve dJ{A)/dA = for A but we can solve it for A: 

4yl5/2 

A = -. ^ J (30) 

3 V^(A - 2)e^/4 erf ( \/A/2j + 2^/A 

The analysis of this expression shows that A(^) exhibits a maximum Ac = 
3.569086042 at Ac = 4.727715383. Therefore there are two variational solutions 
for each < A < Ac, only one for A = Ac and none for A > Ac. This conclusion 
agrees with the rigorous mathematical analysis of the exact solution [5] that we will 
discuss below. Besides, the critical value of the adjustable parameter Ac is also a 
root of d^J{A)/dA^ = 0. 

The exact solution to the one-dimensional Bratu equation (|26p is well-known. 
Curiously enough. He [1], Deeba et al [6], and Khury [7] showed a wrong expression. 
A correct one is (notice that one can write it in different ways) 



.(x)--21n|H^^^ilfclM\ 
{ cosh(0/2) / 



(31) 
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where ^ is a solution to 

^ = cosh(^/2)2 

The critical A-value is the maximum of A(^), and we easily obtain it from the root 
of d\{9)/d9 = that is given by 

e^=(0e-2)-0c-2 = O (33) 

The exact critical parameters are Oo = 2.399357280 and Ac = 3.513830719 that lead 
to u'(0)c = 4. We appreciate that the variational approach provides a reasonable 
qualitative (or even semi quantitative) description of the problem. 

We may try a perturbation approach to the Bratu equation in the form of a 
Taylor series in the parameter A: 

oo 

u{x) =Y,^,{x)\^ (34) 

where, obviously, uo{x) = 0. From the exact expression we obtain 

A A^ A^ 

zi'(O) = + + +...« 0.5A + 0.0417A^ + 0.00625A^ + . . . (35) 
2 24 160 

while the variational approach also yields a reasonable result 
A A^ 43A^ 

u'(0) = - + — + — — + . . . w 0.5A + 0.05A2 + 0.00768A^ + . . . (36) 
2 20 5600 

It seems that the Ritz method already produces satisfactory results for this kind of 

two-point boundary value problems. 

Another simple variational function that satisfies the same boundary conditions 

is 

u{x) = Asm{'Kx) (37) 
It leads to the following variational integral: 

4 2 2 

J{A) = ^-\[I,{A) + L,{A)] (38) 

where Iv{z) and Li,{z) stand for the modified Bessel and Struve functions [8], re- 
spectively. From the minimum condition we obtain 

Att^ 

^ = 2{2 + 7r[/i(A)+Li(A)]} ^^^^ 
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It is not difficult to show that this trial function yields better critical parameters: 
Ac = 3.509329130 and u'(0)c — 3.756549365. Besides, one can easily derive the 
approximate perturbation expansion exactly 



Notice that although the coefficient of A is not exact the remaining ones are more 
accurate than those of the preceding trial function. 

Fig. [T] shows the exact slope at origin u'(0) in terms of A and the corresponding 
estimates given by the two variational functions. We appreciate that the variational 
approach proposed by He [1] yields the solution with smaller u'{0) (lower branch) 
more accurately than the other one (upper branch). This comparison between the 
exact and approximate solutions for a wide range of values of A was not carried 
out before; He [1] simply compared the two slopes at origin for just A = 1. The 
other trial function p7p yields a better overall approximation at the expense of the 
accuracy for small values of A. Some more elaborated approaches, like the Adomian 
decomposition method, fail to provide the upper branch [6]; therefore, the Ritz 
method seems to be suitable for the analysis of this kind of nonlinear problems. 

We may conclude that the Ritz variational method provides a useful insight into 
the Bratu equation. However, one should not forget that there is a relatively simple 
exact solution to this problem and that the generalization of the approach in the 
form of a power series proposed by He [1] 




+ 



(40) 



u{x) = Ax{l — x)(l + CiX + C2X^ + . . .) 



(41) 



may surely lead to rather analytically intractable equations. 



7. Conclusions 

Historically, scientists have developed perturbational, variational and numerical ap- 
proaches to solve nontrivial mathematical problems in applied mathematics and 
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theoretical physics. In some cases, where the exact solution exists but is given 
by complicated special functions, an approximate simpler analytical solution may 
nonetheless be of practical utility. However, He [1] chose examples where either the 
Ritz method does not provide any useful insight, or the exact analytical solutions 
are as simple as the approximate ones, or the direct derivation of the exact result is 
more straightforward than the use of the variational method. From the discussions 
in the preceding sections we may conclude that He's application of the Ritz varia- 
tional method [4] does not show that the approach is suitable for the treatment of 
nonlinear problems. In most of the cases studied here the straightforward analysis 
of the problem yields either more information or the same result in a more direct 
way. 

To be fair we should mention that the Ritz variational method provides a rea- 
sonable bifurcation diagram by means of relatively simple trial functions as shown 
in Fig. [TJ However, even in this case the utility of the approach is doubtful because 
there exists a remarkably simple analytical solution to that equation. The treatment 
of a nontrivial example is necessary to assert the validity of the approach. 

He's choice of the rate equation for chemical reactions [1] is by no means a 
happy one (without mentioning the mistakes in the calculations). In this case the 
exact solution is quite simple and the variational ansatz is unsuitable for practi- 
cal applications. We may argue that a trial function with the correct asymptotic 
behaviour would yield meaningful results. In fact, it may even produce the exact 
result; but one should not forget that such a success would obviously be due to the 
fact that there exist a remarkably simple exact solution available by straightforward 
integration. 

As said before, present results show that He [1] failed to prove that the Ritz vari- 
ational method provides a successful way of treating strongly nonlinear problems. 
Of course, the main ideas behind that variational method are correct, and the case 
of Bratu equation suggests that it may be possible to find appropriate trial func- 
tions for the successful treatment of some problems. Unfortunately, the remaining 
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He's choices [1] do not do much to convince one that the method is worthwhile. 

He's article [1] is an example of the kind of poor research papers that have been 
lately published by some supposedly respectable journals. It is part of such journals' 
policy to reject comments that can reveal this unfortunate situation. One may ask 
oneself what is the profit that those journals get from such a practice. If the reader 
is interested in other examples of poor research papers in supposedly respectable 
journals I suggest some earlier reports in this forum [9-14]. 
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Fig. 1. Bifurcation diagram for the slope at origin u'(0) in terms of A 



